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Abstract 

By identifying non-local effects in systems of identical Bosonic qubits 
through correlations of their commuting observables, we show that 
entanglement is not necessary to violate certain squeezing inequalities 
that hold for distinguishable qubits and that spin squeezing may not 
be necessary to achieve sub-shot noise accuracies in ultra-cold atom 
interferometry. 

1 Introduction 

A generic state of A distinguishable qubits is denned to be entangled if it is 
not fully separable pQ, namely if it cannot be written as 



E(l) „ (2) „ (AO 
PkPk ®Pk • • • ® Pfc 



where the pk > are weights, YlkPk = 1> while each pj^ is a density matrix 
for the j-th qubit acting on the corresponding Hilbert space Hj. 

Entangled A-qubit states have been proposed as means to beat the so- 
called shot-noise limit in metrological applications [2j |3] and spin-squeezing 
techniques have been devised to generate them in systems of ultra-cold 
atoms jH [5] . These states are then used in interferometric experiments where 
their states are rotated by means of collective spin components; in the case 
of distinguishable qubits, all these rotations are local; therefore, preliminary 
spin squeezing of separable states is necessary in order to introduce (non- 
local) quantum correlations. 
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However, in the case of ultra-cold trapped atoms, the qubits involved are 
identical, a fact that asks for a rethinking of the properties and behaviors 
valid for distinguishable qubits. Indeed, the notion of separability based 
on (pQ) is strictly associated with the tensor product structure of the Hilbert 
space, Hat = (&f =1 Hj, which is natural for N distinguishable particles. On 
the other hand, pure Bosonic states must be symmetric under exchange of 
particles and mixed states must be convex combinations of projections onto 
such states. This fact demands a different approach to the notions of non- 
locality and entanglement based not on a structure related to the particle 
aspect of first quantization, as in (JTJ, rather on the behavior of correlation 
functions of commuting observables 19] , more generally related to the 
mode description typical of second quantization. 

In the following, based on a generalized notion of entanglement [S] which 
reduces to the standard one for distinguishable qubits, we show that nei- 
ther entangled states, nor spin-squeezing are necessary in order to achieve 
sub-shot noise accuracies. Though necessary, non-locality comes not from 
the states, rather from the rotations that are implemented in the interfero- 
metric experiments. Indeed, while all rotations are local for distinguishable 
qubits, certain collective spin observables become instead non-local for iden- 
tical qubits and might permit sub-shot noise accuracies without the need of 
a preliminary squeezing of the input state as in [HE]. 



2 Entanglement for Identical Particles 

In general, density matrices as in ([T]) are not allowed Bosonic states not even 
if p^ 1 = pk for all i. Indeed, consider two qubits and fix an orthonormal 
basis C 4 3 i,j =t, t- If they are Bosons, their states cannot have 

I It) - I T4) 

non- vanishing components on the anti-symmetric state |^_) = —j= . 

Since ($-\p <8> p\^?-) = Det(p), the density matrix p <g> p cannot correspond 
to a state of two Bosonic qubits unless p is a projection (Det(p) = 0). 

Therefore, the tensor product structure which is natural for distinguish- 
able particles is not appropriate for discussing the entanglement properties 
of systems of indistinguishable particles. These should rather be investigated 
within the second quantization formalism whereby one introduces creation 
and annihilation operators af of single-particle orthonormal basis states |i), 
obeying the canonical commutation relations [a, , aj] = 5%j. 

Entanglement in such a context should correspond to whether, given a 
state to of the system, there are non-classical correlations among commuting 
observables 012]; for instance, between non-overlapping spatial regions Vi 2 . 
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Specifically, this can be inspected by considering the structure of two-point 
functions of the form u(Pv 1 Pv 2 )j where Py a is any polynomial in creation 
and annihilation operators a^(ip a ) of states \ip a ) spatially localized in the 
region V a . The relevant fact is that the spatially local algebras generated by 
polynomials Py a commute: [a(ipi) , 0^(1^2)] = (V'ilV'2) = 0. 

More in general, one can argue about the entanglement between observ- 
ables belonging to two generic commuting sub-algebras (A, B) of the entire 
algebra generated by creation and annihilation operators, which we shall re- 
fer to as algebraic bipartition [§J. We shall call an operator (A, S)-local if 
of the form AB, A 6 A and B 6 B and a state u (A, £>)-separable if the 
expectations u(AB) of local operators AB can be decomposed into a convex 
linear combination of product of expectations: 

u(AB) = X ^t(A) cu^B) , A, > , ]T A, = 1 , (2) 

i i 

in terms of other states u^ ,b ; otherwise, u is (A, i3)-entangled. 

Remark 1. In the case of two qubits, the above definitions reproduce the 
standard notions if one chooses the algebraic bipartition A = B = M2, 
where M2 the algebra of 2 x 2 matrices over C 2 , and u(AB) = Tr (p A cg> B), 
with p a two-qubit density matrix. However, for identical particles, there 
is no a priori given bi-partition so that questions about entanglement and 
separability, non-locality and locality are meaningful only with reference to 
a specific class of (commuting) observables. 



2.1 TV Bosons in a Double- Well Potential 

A concrete application of the previous considerations is the second quantiza- 
tion of a single-particle with Hilbert space C 2 which, in the Bose-Hubbard 
approximation, effectively describes iV ultra-cold atoms confined by a double- 
well potential. Then, the state | \) describes one atom located within the 
left well and the state | f) an atom localized within the right one. Let |0) be 
the vacuum state and a\ W the creation operators of a particle in the states 
I 4) and I |), that is a*\0) = | |), &t|0) = | t). 

When the total number iV is conserved, the symmetric Fock space of 
this two-mode system is generated by iV + 1 orthonormal eigenvectors of the 
number operator a^a + tfb: 



a, 1 



'^7ira |0> ' °- k - N - (3) 
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Because of the orthogonality of the spatial modes, by considering the norm- 
closures of all polynomials P a in a,a^, respectively P& in b,b', one obtains 
two commuting subalgebras A and B. 

According to the states \k) are (A, B) -separable; indeed, they are 
created by the (A, £>)-local operators (a^) k (b j( ) N ~ k . More in general, (A,B)- 
separable states must be convex combinations of projections \k){k\ [S]: 

N N 

P = ^Pk\k)(k\ p k >0,J2Pk = l- (4) 

k=0 k=0 

Consider instead the following operators 

J x = -(a f & + ab ] ) , J y = ^(a f 6 - atf) , J z = ^(a f a - tfb) , (5) 



that satisfy the SU(2) algebraic relations [J x , J y ] = i J z . They are all non- 
local with respect to the algebraic bipartition (A, B) and such are the expo- 
nentials e WJ * and e iflJ » , while e i9Jz = e i0a+a e^ 6t6 i s S)-local. 

By means of a Bogolubov transformation to other creation and annihi- 
lation operators (c # ,<i # ), such that a = and b = — — , one obtains 
another bipartition (C, T>) and rewrites 

Jx = - (c^c — d)d) , J y = ~^~.{d)c — cfct) , J z = - (c^d + cd)) . 

Relatively to {C,V}, it is now e i9Jx = e ie ^%~^ d ^ d w hi c h acts locally 

Remark 2. In first quantization, an (^4., inseparable state for a iV = 2 

I +| \ _|_ I ||\ 

Bosonic qubits like |11) = a^b^O) corresponds to -j= . Such state 

is surely entangled for distinguishable qubits, while, according to our defi- 
nition, it is no longer so for identical Bosonic qubits; the reason is that its 
entanglement is only formal as it comes from the necessary symmetrization 
of the separable state | tl) [ 10J . 

A Bogolubov transformation as the one above corresponds to a change of 
basis in the single particle Hilbert space, from the one of spatially localized 

states, to the one of c + |0) = -=(\ |) + | |», d f |0) = -j=(\ |) - | t})- Physi- 

v 2 v 2 

cally speaking, such states are eigenstates of the single particle Hamiltonian 

in the Bose-Hubbard approximation with a highly penetrable barrier. The 

change to the energy bipartition (C, T>) is non-local with respect to the spatial 

bipartition (A,B), though it corresponds to a local unitary rotation in first 

quantization. 
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2.2 Collective Spin Inequalities and Entanglement 

In the case of a system of N distinguishable qubits, the collective angular 
momentum operators J x ,y,z and the corresponding rotations are sums of single 

particle spin operators, Jxjy,*, i- e. J x , y ,z = ^2f=i Jx,y,z- These operators are 
local with respect to the tensor product structure in (JJJ). 

Based on this, the variance A 2 of the collective spin = n x J x + 
n y J y + n z J z along the unit spatial direction n = (n x , n y , n z ), with respect to 
separable vector states |^) = 0^=1 \4>j) results 



A 2 J H = (*| J*\V) - (*| J B |tt> 5 

N 

A £(toi4%>) a < T - ( 6 ) 

Therefore, A 2 is an entanglement witness for pure states, in the sense that 
if A 2 Jrf > N/4 then the pure state cannot be fully separable. 

This is no longer the case for N identical Bosonic qubits. Indeed, consider 
the number states \k) in (jHj); using ([5]), one gets 

(k\J H \k) = ^(2k-N) (7) 



N + 2k(N-k) 



4 

2 N(N - 1) -Gk(N- k) 

Hz 4 
1-n 2 



4\k) 



A 2 J H = [N + 2k(N - k)j . (9) 

Therefore, if k ^ 0, N, for all n that satisfy 

n 2 < 2k{N ~ k) < 1 
z N + 2k(N - k) ~ ' 

the states though (^4, £>)-separable, nevertheless yield A 2 .1^ > N/4; 
therefore, A 2 is not an entanglement witness for pure states of Bosonic 
qubits. 

In greater generality, inequalities for mean values and variances of col- 
lective spin operators with respect to any (mixed) separable state of distin- 
guishable qubits ([T]) have been derived in [FT]; these are called spin squeezing 
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inequalities and read 



<4> + <4> + <4> - < o . (io) 

A 2 J Hl + A 2 + A 2 J Hs -j>0, (11) 
(4) + (4) - y - ( N - ^Jn, < , (12) 



(A- - l)(A 2 J, 7j - A-.7, 7 J - (4) - N{N 2) > . (13) 



where ni,2,3 denotes any triplet of unit vectors corresponding to orthogonal 
spatial directions and (X) the mean value of an operator X. 

It is thus interesting to study whether these inequalities are also satisfied 
by (A, £>)-separable states of N identical qubits. Let (k a ) = Y^k=oPk^ a -> 
a = 1,2, denote first and second moments of the N + 1- valued stochastic 
variable k with respect to the probability distribution n = {pk}k=o- Us- 
ing ([^), (|SJ), mean- values and variances of collective spin operators J a with 
respect to the states in (jlj) read 

(Jn) = %(2(k)-N) (14) 

2\ 



\Jfi) 



A 2 X n 



2 

N{l + 2(k))-2(k 2 ) 
4 

,2 



+ ^{N{N -l)-QN{k) + Q{k 2 )^ (15) 
N{1 +2(k)) ~2(k 2 ) 



4 

,2 



+ ^(6(A; 2 )-2(fc)(iV + 2(A;))-iv) . (16) 

From the orthogonality of the triplet r? lj2 ,3, it follows that ^ 2 2 + n L + n iz = 1; 
one can thus check that all inequalities but f lT2"j) are satisfied by (A, B)- 
separable states. Concerning ( jUJ) , its left hand side reads 

A/ A 2 k- (k)(N- (k)) 



+ ^fi K (N + 2)(k)(N-(k))-3NA 2 k ), (17) 

where A 2 k := (k 2 ) — (k) 2 is the variance of k with respect to 7r = {paJ^Lo- If 7r 
is chosen such that a := (k)(N — (k)) > NA 2 k, then 5 becomes positive and 



1 In [11) . these inequalities are derived with respect to the standard triplet rC\ = x, 
H2 — y, «3 = z. The result easily extends to more general triplets. 
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thus f[T2l results violated by the corresponding (A, inseparable states (TjJ 
for all orthogonal triplets with 



1 > < > 



N(a - A 2 k) 



(18) 



(N + 2)a-3N A 2 k ' 



Consider the pure states \£) in with £ ^ 0,N; in such a case, pk = 6^, 
(k) = £ and A 2 k = 0, so that (JI2D is violated for 1 > n 2 z > N/(N + 2). 

3 Spin Squeezing and Metrology 

The preceding results indicate that spin-squeezing inequalities that are de- 
rived for distinguishable qubits can not directly be used as entanglement 
witnesses in the context of identical qubits. Since the use of spin-squeezed 
states for metrological purposes have recently become the focus of much the- 
oretical [21 1121 H3] investigations, we now discuss the impact of particle 
indistinguishability on such an issue. 

3.1 Spin Squeezing 

For any orthogonal triplet of space-directions ^1,2,3, the Heisenberg uncer- 
tainty relations for the SU(2) operators read 



One speaks of spin-squeezing when one of the variances can be made smaller 



otherwise unavailable accuracies has been studied in relation to the measure 
of an angle 9 by interferometric techniques. These are based on a rotation of 
an input state p into 



and upon measuring on pg the collective spin J^ 2 , where n,2 J- fix. By choosing 
the remaining orthogonal unit vector n 3 such that (Ja 3 ) = Tr(pJfj 3 ) 7^ 0, 
by error propagation, the uncertainty 58 in the determination of 9 can be 
estimated by [3] 



(19) 



than i 



(Jn 3 ) ■ The relevance of states satisfying this condition for achieving 




(20) 



5 2 9 




(21) 
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in terms of the spin-squeezing parameter 



9 NA 2 J H , , , 

* : ^ W ' (22) 

The value = — is called shot-noise limit; in the case of distinguishable 

N 

qubits, it gives the lower bound to the attainable accuracies when the input 
state p is separable. Indeed, in such a case one finds < 1. This result 
follows from the inequality £yy — €i where the new spin-squeezing parameter 

e a ;= N&J* ( (23) 

has been introduced in [12]; by means of the local structure of the collective 
spin operators J x , y ,z, one can prove that £f is always > 1 for separable states 
of distinguishable qubits. Therefore, using distinguishable qubits, the shot- 
noise limit can be beaten, namely accuracies better than 1/N can be achieved 
only if £| < 1, that is only by means of entangled states. 

Let us instead consider N identical Bosonic qubits in the (.4, i3)-separable 
pure states \k) and any triplet of orthogonal spatial directions ^1,2,3 with 
ni 7^ z. Using ©-© an d n\ z + n| a + n\ z = 1, one computes 



NA 2 J Hl AT l-n\ z N + 2k(N-k) 



^ (Jn 2 y + (Jn 3 y nl + nl (2k - N) 2 

N(N + 2k(N — k)) , , . 

= ( 2k -Nr - 1 • °- k - N - (24) 

In the case of the (A, inseparable density matrices (j3j^, first observe that, 
thanks to the Cauchy-Schwartz inequality, one has J2k=o Pk(Jn)k > (Jn) 2 , 
where (X)\ denotes the mean-value of X with respect to the number state 
\k). Then, 

£2 y J2k=0 Pk ^fc^ni 

2 EtoPkN( y N + 2k(N-k) 



> 



U2z + n3z (EtoPk(N-2k)) 

Ek=oPk (fc- f) 2 
(EL,** (^-f) N2 ' 



2 
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where the last inequality follow from the second line of (|24l) . A further 
application of the Cauchy-Schwartz inequality to the right hand side of the 
last inequality yields £§ > 1 for all (A, ^-separable states when rii ^ z. 

If one chooses ni = z, in the case of (A, ^-separable mixed states, one 
finds A 2 J z 7^ and (Jft 2 ,n 3 ) — 0; therefore, £| (and £j^) diverges. Instead, for 
(A, inseparable pure states |fc) also A k J z = whence £| (and £j^) are not 
defined and must thus be computed by means of suitable limiting procedures. 

Let us consider the (A, i3)-entangled vector state |^) = J2k=o \fv~k \k) , 
with real coefficients from a probability distribution n = {pk\k=o over the 
stochastic variable k. Then, from ([I])-© it follows that ( J y ) = and A 2 J z = 
A 2 k; therefore £f = and 

2 _ NA 2 J Z _ NA 2 k 

£w ~ I T \2 ~~ 7 ^ \~2 ■ \ ^) 

{ x) ( £j=i y/k{N -k+l) y/pipi=i) 
In the case of a Gaussian distribution peaked around k = £ ^ 0, N, 

p k = - exp — - , Z = }^exp — - , (26) 



7 fe=o 

one finds 

2iV + O fe~^ 

& = ~ * '- ~ 2 ■ (27) 

(y/(e + i)(N-e) + y/e(N-e + i)) 

Thus, for sufficiently small a, $y < 1 for all £ ^ 0,N. 
On the other hand, by choosing 

AT 

l*> = £ E l*> + (l-P)l^/2>. 0<p<l, (2f 

k=0,k^N/2 

it turns out that A 2 /c = p(N + 2)(N+ 1)/12 so that (T25D yields 

iV(iV + l) 
^ ^(v/T^ + g) 2 ' 



where 
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P 



iV 2 (iV + 2) 

v ' k^N/2,N/2+l 



VKN-k + l) . 



Letting p ->• one gets |-0) ->■ |iV/2) and, if iV > 3, 

e w ^N{N + l)/l2> 1 . 
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Remark 3. The above two examples show that, when n\ = z, the spin- 
squeezing parameters £w,s are not well-defined: different values for £j^ 5 
can be obtained by approaching a state \k) via different limit procedures. 
This fact is also of practical importance: indeed, in |T3j, approximations to 
Fock states \k) have been experimentally constructed that are characterized 
by spin-squeezing parameters < 1. This property arises from the fact 
that the approximations are [A, £>)-entangled states. The previous discus- 
sion shows that some care has to be taken in constructing the perturbations 
of \k); indeed, not all (A, £>)-entangled states arbitrarily close to it automat- 
ically have £yys — 1- Therefore, the spin-squeezing parameters are not 
always useful for metrological applications, a better quantity is the so-called 
quantum Fisher information [T3] , which as we shall show below, is continuous 
and well defined for all Bosonic qubits. 



3.2 Quantum Fisher Information 

In a measurement of the angle 9 based on the state rotation f[2"Uj) . the error 
AO given by a locally unbiased estimator E of the angle 6 is bounded by (see 
the Appendix) 

A 2 6 > — ^— , (29) 

where F[p, J^J is the so-called quantum Fisher information associated with 
the rotation of p around n\. 

In order to overcome the shot-noise limit A 2 9 = 1/N, the quantum Fisher 
information must then be strictly larger than N. In the appendix it is also 
showed that, in full generality, 

F[p,J Hl }A 2 J H2 >(J H3 ) 2 , (30) 

where n.1,2,3 is a triplet of orthogonal spatial directions. Thus, if (Jn 3 ) 7^ 0, 
one gets the following relation between the quantum Fisher information and 
the squeezing parameter £jy in (T22|) : 



F ^^^ = l' (31) 

In the case of distinguishable qubits, from (129]) and (13TT) it follows that spin- 
squeezing, namely < 1, opens the possibility of achieving A 2 9 < 1/N, 
thus of beating the shot-noise limit. 

In the case of identical qubits and of (A, £>)-separable states, the right 
hand side of the above inequality diverges if Hi — z as (Jn 3 ) = 0, while it 
does not make sense if n 2 = z for then also A 2 J^ 2 = whence, as already 
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observed, £jy is not defined. However, the quantum Fisher information is 
always well-defined. Indeed, using f[3"5j) and OH]), one finds that, if k ^ 0,N, 

F[\k)(k\,J ri }=4A 2 J z = (l-n 2 z )^N + 2k(N-k)^ >N (32) 

2 2k(N-k) 
for < n z < N + 2k ^ N _ k) < L In particular, 

F[\k)(k\, J y ) = N + 2k(N-k) , 

so that, according to f[2"9j) . for all k ^ 0, A" the (.4., inseparable pure states \k) 
might overcome the shot-noise limit, with the twin Fock state \N/2) yielding 
F[p, J y ] = 0(N 2 ) thus permitting to approach the so-called Heisenberg limit 
A 2 6 = l/N 2 . 

Remark 4. Notice that even approximating a number state \£) by an ex- 
perimentally more realistic superposition \£, a) of states \k) with coefficients 
as in (|2T)|) . may beat the shot noise limit. Indeed, one computes 

F[\£, a)(£, a\, J y ] = 4 A 2 J y = N + 2£(N - £) + O (e~^ 

which can be kept > iV by suitably small a. 

Instead, making the quantum Fisher information larger than iV is impos- 
sible without (A, B)-non-locality; indeed, F[\£)(£\,J Z ] = AA 2 J Z = 0. Even 
considering the (A, £>)-entangled perturbation \£, a) does not help; indeed, 



F[\£,a)(£,a\,J z ] = AA 2 J Z = 8e"^ + O (e" 



Therefore, the lower bound to the error in ( |29|) becomes arbitrarily large 
when a — > 0. 

When dealing with (A, £>)-separable mixed states (jl]), by means of equa- 
tion ( 13T|) in the Appendix, one computes [H] 

F\p,J H ] = (l-n 2 z )(N + 2N(k)-(k 2 ) (33) 

N 

- *E^<*+i><"-*>)- w 

Thus, if ( |32|) holds for a certain \£), then, by continuity, F[p, Jg] > A^ for a 
probability distribution n = {pk)k=o suitably peaked around k = £, hence 
able to overcome the shot-noise limit. 

On the other hand, from the previous section we know that for all such 
mixed states > £| > 1; therefore, based on this lower bound to the 
squeezing parameter, we would wrongly discard such states as not useful for 
metrological applications. 
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4 Discussion 



Unlike for distinguishable qubits, in the case of identical Bosonic qubits, en- 
tangled states are not necessary to reach sub-shot noise accuracies in param- 
eter estimation. This phenomenon is surely due to the non-local character 
of the system; however, the non-locality is not in the states \k) which in fact 
are (A, inseparable, rather in the (A, £>)-non-local character of the state- 
rotation ( 120]) generated by J ft ^ J z . Were the qubits distinguishable, neither 
the state nor the rotation would carry elements of non-locality so that in or- 
der to beat the shot-noise limit, the state should be turned into an entangled 
one before feeding the interferometric apparatus practically implementing 
the state- rotation (1201) . This is exactly what is done via a spin-squeezing 
technique in the experiments reported in [H |5]. Instead, the main point we 
make here is that in experiments involving identical qubits, no preliminary 
squeezing is needed before rotating the state. One might as well do with, 
say, a state |iV/2) , the rotation around n ^ z taking care of introducing the 
necessary non-locality. 



5 Appendix 



A most used quantum Fisher information F[p, Jo] is given by [16] 

F[p,J n ] :=tr(pL 2 ) , (35) 



where L, known as symmetric logarithmic derivative, is a Hermitean operator 
such that 

__ r P -±^ = -i[Jn,p}. (36) 
Given a spectral decomposition p = J2j r j l r j)( r jl' one computes 



d, 



F[p,J H } = 2 



Ti - Ta 



yi \Jfi yj) 



(37) 



From such an expression one sees that the quantum Fisher information is a 
continuous function of the state p and that, for pure states, 



Jfi, 



(38) 



1; then, inequality 



An estimator E is locally unbiased if dgTripo E) 

6»=u 

follows from applying to this relation the Cauchy-Schwartz inequality for 
matrices 



Tr(AB) < Tr(A f A)Tr(5 t B) 
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Analogously, inequality ( 15D|) follows from the fact that Tr(p e ( Jf l2 — (Jn 2 )e) = 
implies 

References 

[1] R. Horodecki et al., Rev. Mod. Phys. 81, 865 (2009) 

[2] M. Kitagawa, M. Ueda, Phys. Rev. A 47, 5138 (1993) 

[3] D.J. Wineland, J. J. Bollinger, W.M. Itano, D.L. Heinzen, Phys. Rev. 
A 50, 67 (1994) 

[4] C. Gross et al, Nature 464, 1165 (2010) 

[5] M.F. Riedel et al, Nature 464, 1170 (2010) 

[6] P. Zanardi, Phys. Rev. A 65, 042101 (2002) 

[7] H. Narnhofer, Phys. Lett. A 310, 423 (2004) 

[8] L. Viola et al, Contemp. Math. 38, 117 (2205) 

[9] F. Benatti, R. Floreanini, U. Marzolino, Ann. Phys. 325, 924 (2010) 

[10] G. Ghirardi, L. Marinatto, T. Weber, J. Stat. Phys. 108, 49 (2002) 

[11] G. Toth et al., Phys. Rev. A 79, 042334 (2009) 

[12] A. S0rensen et al, Nature 409, 63 (2001). 

[13] L. Pezze, A. Smerzi, Phys. Rev. Lett. 102, 100401 (2009) 

[14] J. Esteve et al, Nature 455, 1216 (2008) 

[15] S.L. Braunstein, CM. Caves, G.J. Milburn, Annals of Physics 247, 135 
(1996) 

[16] S. Luo, Lett. Math. Phys. 53, 243 (2000) 



13 



